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Abstract: This paper provides necessary and sufficient conditions for constructing a uni- 
versal quantum computer over continuous variables. As an example, it is shown how a 
universal quantum computer for the amplitudes of the electromagnetic field might be con- 
structed using simple linear devices such as beam-splitters and phase shifters, together 
with squeezers and nonlinear devices such as Kerr-effect fibers and atoms in optical cav- 
ities. Such a device could in principle perform 'quantum floating point' computations. 
Problems of noise, finite precision, and error correction are discussed. 

Quantum computation has traditionally concerned itself with the manipulation of 
discrete systems such as quantum bits, or 'qubits'. 1-2 Many quantum variables, such 
as position and momentum, or the amplitudes of electromagnetic fields, are continuous. 
Although noise and finite precision make precise manipulations of continuous variables 
intrinsically more difficult than the manipulation of discrete variables, because of the recent 
developments in quantum error correction 3-5 and quantum teleportation 6-7 of continuous 
quantum variables it is worthwhile addressing the question of quantum computation over 
continuous variables. 

At first it might seem that quantum computation over continuous variables is an 
ill-defined concept. First consider quantum computation over discrete variables. A uni- 
versal quantum computer over discrete variables such as qubits can be defined to be a 
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device that can by local operations perform any desired unitary transformation over those 
variables. 1-2 ' 8 More precisely, a universal quantum computer applies 'local' operations 
that effect only a few variables at a time (such operations are called quantum logic gates): 
by repeated application of such local operations it can effect any unitary transformation 
over a finite number of those variables to any desired degree of precision. Now consider the 
continuous case. Since an arbitrary unitary transformation over even a single continuous 
variable requires an infinite number of parameters to define, it typically cannot be approx- 
imated by any finite number of continuous quantum operations such as, for example, the 
application of beam-splitters, phase-shifters, squeezers, and nonlinear devices to modes of 
the electromagnetic field. It is possible, however, to define a notion of universal quantum 
computation over continuous variables for various subclasses of transformations, such as 
those that correspond to Hamiltonians that are polynomial functions of the operators cor- 
responding to the continuous variables: a set of continuous quantum operations will be 
termed universal for a particular set of transformations if one can by a finite number of 
applications of the operations approach arbitrarily closely to any transformation in the 
set. 

This paper provides necessary and sufficient conditions for universal quantum compu- 
tation over continuous variables for transformations that are polynomial in those variables. 
Such a continuous quantum computer is shown to be capable in principle of performing 
arithmetical manipulations of continuous variables in a 'quantum floating point' compu- 
tation. Issues of noise and finite precision are discussed and applications proposed. 

Consider a single continuous variable corresponding to an operator X. Let P be the 
conjugate variable: [X, P] = i. For example, X and P could correspond to quadrature 
amplitudes of a mode of the electromagnetic field (the quadrature amplitudes are the 
real and imaginary parts of the complex electric field). First investigate the problem of 
constructing Hamiltonians that correspond to arbitrary polynomials of X and P. It is 
clearly necessary that one be able to apply the Hamiltonians ±X and ±P themselves. In 
the Heisenberg picture, applying a Hamiltonian H gives a time evolution for operators 
A = i[H, A], so that A(t) = e iHt A(0)e~ iHt . Accordingly, applying the Hamiltonian X for 
time t takes X — > X, P — > P — t, and applying P for time t takes X — > X + 1, P — > P: the 
Hamiltonians X and P have the effect of shifting the conjugate variable by a constant. In 
the case of the electromagnetic field, these Hamiltonians correspond to linear displacements 
or translations of the quadrature amplitudes. 

To construct an arbitrary Hamiltonian of the form aX + bP + c, first apply X for a 
time adt, where dt is a short period of time, then apply P for a time bdt, and finally apply 
P for a time V cdt, X for a time V cdt, —P for a time ^fcdt, and —X for a time V cdt. The 
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net effect is a transformation 

g— iXV cdt^—iPVcdtgiXV cdt^iPV cdt ^ibPdt ^iaXdt 

= 1 + i( - i[X, P]c + bP + aX)dt + 0(dt 3/2 ) (1) 

^ A(aX+bP+c)dt 

By making dt sufficiently small, one can approach arbitrarily closely to effecting a Hamil- 
tonian of the desired form over small times. By repeating the small-time construction t/dt 
times, one can approach arbitrarily closely to effecting the desired Hamiltonian over time 
t. 

There are clearly simpler ways to enact an overall phase shift. When applied to 
arbitrary sets of Hamiltonians, however, the construction given above provides the pre- 
scription for determining exactly what Hamiltonians can be constructed by the repeated 
application of Hamiltonians from the set: if one can apply a set of Hamiltonians {±Pj}, 
one can construct any Hamiltonian that is a linear combination of Hamiltonians of the 
form ±i[Hi, Hj], ±[H i7 [Hj,H k ]}, etc., 9-13 and no other Hamiltonians. That is, one can 
construct the Hamiltonians in the algebra generated from the original set by commuta- 
tion. This key point, originally derived in the context of quantum control and discrete 
quantum logic, makes it relatively simple to determine the set of Hamiltonians that can 
be constructed from simpler operations. 

The application of the translations ±X and ±P for short periods of time clearly allows 
the construction of any Hamiltonian that is linear in X and P; this is all that it allows. 
Suppose now that one can apply the quadratic Hamiltonian H = (X 2 + P 2 )/2. Since 
P = i[H,P] = X, X = i[H,X] = — P, application of this Hamiltonian for time t takes 
X — > costX — siniP, P — > costP + sintX. If X and P are quadrature amplitudes of 
a mode of the electromagnetic field, then H is just the Hamiltonian of the mode (with 
frequency u = 1) and corresponds to a phase shifter. Hamiltonians of this form can be 
enacted by letting the system evolve on its own or by inserting artificial phase delays. Note 
that since e lHt is periodic with period 1/4-7T, one can effectively apply —H for a time St 
by applying H for a time 4rr — St. The simple commutation relations between P", X and 
P imply that the addition of ±H to the set of operations that can be applied allows the 
construction of Hamiltonians of the form aH + bX + cP + d. 

Suppose that in addition to translations and phase shifts one can apply the quadratic 
Hamiltonian ±S = ±(XP + PX)/2. S has the effect X = i[S,X] =X,P = i[S, P] = -P, 
i.e., applying +S takes X — * e t X, P — > e~*P: S 'stretches' X and 'squeezes' P by 
some amount. Similarly — S squeezes X and stretches P. In the terminology of quantum 
optics, S corresponds to a squeezer operating in the linear regime. It can easily be verified 
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that [H, S] = i(X 2 — P 2 ). Looking at the algebra generated from X,P,H and S by 
commutation, one sees that translations, phase shifts, and squeezers allow the construction 
of any Hamiltonian that is quadratic in X and P, and of no Hamiltonian of higher order. 

To construct higher order Hamiltonians, nonlinear operations are required. One such 
operation is the 'Kerr' Hamiltonian H 2 = (X 2 + P 2 ) 2 , corresponding to a % 3 process in 
nonlinear optics. This higher order Hamiltonian has the key feature that whereas com- 
muting the previous Hamiltonians, X, P, H, S with some polynomial in X and P resulted 
in a polynomial with the same or lower order, commuting H 2 with a polynomial in X and 
P typically increases its order, e.g., 

[H 2 , X] = H[H, X] + [X, H]H = i(HP + PH) = i(X 2 P + PX 2 + 2P 3 ) /2 (2.1) 
[H 2 ,P] = -i(P 2 X + XP 2 + 2X 3 )/2 (2.2) 
[H 2 ,S] = H[H,S] + [H,S]H = i(X 4 - P 4 )/2 . (2.2) 

By evaluating a few more commutators, e.g., [X, [H 2 , S}] = P 3 , [P, [H 2 , S]] = X 3 one sees 
that the algebra generated by X, P, H, S and H 2 by commutation includes all third order 
polynomials in X and P. A simple inductive proof now shows that one can construct Hamil- 
tonians that are arbitrary Hermitian polynomials in any order of X and P. Suppose that 
one can construct any polynomial of order M or less, where M is of degree at least 3. Then 
since [P 3 ,P m X n ] = iP m+2 X n - 1 + lower order terms, and [X s , P m X n ] = iP m - 1 X n+2 + 
lower order terms, one can by judicious commutation of X s and P 3 with monomials of 
order M construct any monomial of order M + 1. Since any polynomial of order M + 1 can 
be constructed from monomials of order M + 1 and lower, by applying linear operations 
and a single nonlinear operation a finite number of times one can construct polynomials 
of arbitrary order in X and P to any desired degree of accuracy. Comparison with similar 
results for the discrete case 14 shows that the number of operations required grows as a 
small polynomial in the order of the polynomial to be created, the accuracy to which that 
polynomial is to be enacted, and the time over which it is to be applied. 

The use of the Kerr Hamiltonian H 2 was not essential: any higher order Hamiltonian 
will do the trick. Note that commutation of a polynomial in X and P with X and P 
themselves (which have order 1) always reduces the order of the polynomial by at least 
1, commutation with H and S (which have order 2) never increases the order, and com- 
mutation with a polynomial of order 3 or higher typically increases the order by at least 
1. Judicious commutation of X,P,H and S with an applied Hamiltonian of order 3 or 
higher therefore allows the construction of arbitrary Hermitian polynomials of any order 
in X and P. 
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The above set of results shows that simple linear operations, together with a single 
nonlinear operation, allow one to construct arbitrary polynomial Hamiltonian transfor- 
mations of a single quantum variable. Let us now turn to more than one variable, e.g., 
the case of an interferometer in which many modes of the electromagnetic field inter- 
act. Suppose now that there are many variables, {Xi,Pi}, on each of which the simple 
single-variable operations described above can be performed. Suppose in addition Hamil- 
tonians of the form = ±(PiXj — XiPj) can be applied. Since = i[Bij,Xj\ = Xj 

Xj = i[P>ij, Xj] = —Xi, Pi = i[Bij,Pi] = Pj, Pj = i[Bij, Pj] = —Pi, this operation has the 
effect of taking Ai — > cos tA t + sin tAj , Aj — > cos tAj — sin tAi , for A t = P t , Aj = Xj , Pj . 
For the electromagnetic field, Bij functions as a beam splitter, linearly mixing together the 
two modes % and j. By repeatedly taking commutators of Bij with polynomials in Xi, Pi, 
it can be easily seen by the same algebraic arguments as above that it is possible to build 
up arbitrary Hermitian polynomials in {Xi, Pi}. 

This concludes the derivation of the main result: simple linear operations such as 
translations, phase shifts, squeezers, and beam splitters, combined with some nonlinear 
operation such as a Kerr nonlinearity, suffice to enact to an arbitrary degree of accuracy 
Hamiltonian operators that are arbitrary polynomials over a set of continuous variables. 
Note that in contrast to the case of qubits, in which a nonlinear coupling between qubits 
is required to perform universal quantum computation, in the continuous case only single 
variable nonlinearities are required, along with linear couplings between the variables. 

In analog with information over classical continuous variables, which is measured in 
units of 'nats' (1 nat = log 2 e bits), the unit of continuous quantum information will be 
called the 'qunat.' Two continuous variables in the pure state | -0) 12 possess — trpilnpx 
qunats of entanglement, where p\ = tr 2 IV 7 ) 12 I • F° r two squeezed vacua (squeezed by an 
amount e~ r ) entangled using a beam splitter as in refs. (5-7) the entropy so computed 
from the approximate EPR state is given by 

S(p) = (1 + n) ln(l + n) — n\nn qunats (3) 

with n = e r sinhr. For example, e 2r = 10 gives 10 dB of squeezing in power, corresponding 
to r = 1.15129. By equation 3, two continuous variables entangled using a 10 dB squeezer 
then possess 2.60777 qunats of shared, continuous quantum information, equivalent to 
3.76221 qubits of discrete quantum information. 

Quantum computation over continuous variables can be thought of as the systematic 
creation and manipulation of qunats. Universal quantum computation for polynomial 
transformations of continuous variables effectively allows one to perform 'quantum floating 
point' manipulations on those variables. For example, it is clearly possible using linear 
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operations alone to take the inputs Xi, X 2 and to map them to Xi, aXi+bX 2 +c. Similarly, 
application of the three- variable Hamiltonian X1X2P3 allows one to multiply X\ and X 2 
and place the result in the 'register' X 3 : 

X t = i[X 1 X 2 P 3 ,X 1 ] = 0,X 2 = iiX^P^Xs] = 0,X 3 = {[XiXtP^Xz] = x x x 2 

(4) 

X\ Xi, X 2 — > X 2 , X% — > X% + X\X 2 t. 

A wide variety of quantum floating point operations are possible. Any polynomial trans- 
formation of the continuous variables is clearly possible, as is any transformation that 
can be infinitesimally represented by a convergent power series. Just as classical compu- 
tation over continuous variables in principle allows one to solve problems more rapidly 
than is possible digitally, 14 it is interesting to speculate that quantum computation over 
continuous variables might in principle allow the solution of problems more rapidly than is 
possible using a 'conventional,' discrete quantum computer. Continuous variable compu- 
tation has its own set of problems that might be sped up by the application of continuous 
quantum computation: for example, such a continuous quantum computer might be used 
to investigate continuous A^P-complete problems such as the 4-Feasibility problem, that 
is, the problem of deciding whether or not a real degree 4 polynomial in n variables has a 
zero. 15 In practice, of course, due to finite precision a continuous quantum computer will 
effectively be able to solve the same set of problems that a 'conventional' discrete quantum 
computer can, although it may be able to perform some operations more efficiently. 

The ability to create and manipulate qunats depends crucially on the strength of 
squeezing and of the nonlinearities that one can apply. 10 dB squeezers (6 dB after at- 
tenuation in the measurement apparatus) currently exist. 16 High Q cavity quantum elec- 
trodynamics can supply a strong Kerr effect in a relatively lossless context, and quantum 
logic gates constructed for qubits could be used to provide the nonlinearity for continuous 
quantum variables as well. 17 Here the fact that only single- mode nonlinearities are required 
for universal quantum computation simplifies the problem of effecting continuous quantum 
logic. Nonetheless, the difficulty of performing repeated nonlinear operations in a coherent 
and loss-free manner is likely to limit the possibilities for quantum computation over the 
amplitudes of the electromagnetic field. 

Noise poses a difficult problem for quantum computation, 18-20 and continuous vari- 
ables are more susceptible to noise than discrete variables. Since an uncountably infinite 
number of things can go wrong with a continuous variable, it might at first seem that 
continuous error correction routines would require infinite redundancy. In fact, continuous 
quantum error correction routines exist and require no greater redundancy than conven- 
tional routines. 3-5 Such routines are capable of correcting for noise and decoherence in 
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principle: in practice, measurement noise, losses, and the lack of perfect squeezing will lead 
to imperfect error correction. 5 Surprisingly, continuous quantum error correction routines 
are in some sense easier to enact than discrete quantum error correction routines, in that 
the continuous routines can be implemented using only linear operations together with clas- 
sical feedback. 5 The relative simplicity of such routines suggests that robust, fault-tolerant 
quantum computation may in principle be possible for continuous quantum variables as 
well as for qubits (A scheme for quantum computation is fault-tolerant if quantum com- 
putations can be carried out even in the presence of noise and errors. 21-22 A fault-tolerant 
scheme that allows for arbitrarily long quantum computations to be carried out is said 
to be robust. 23 ). If this is indeed the case then quantum computation over continuous 
variables, despite its intrinsic difficulties, may be an experimentally viable form of quan- 
tum information processing. Continuous variables might be used to simulate continuous 
quantum systems such as quantum field theories. Even in the absence of fault tolerance, 
the large bandwidths available to continuous quantum computation make it potentially 
useful for quantum communications and cryptography. 24 
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